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Charged particles (X) decaying after primordial nucleosynthesis are constrained by the require- 
ment that their decay products should not change the light element abundances drastically. If the 
decaying particle is negatively charged (X~ ) then it will bind to the nuclei. We consider the effects 
of the decay of X when bound to Helium-4 and show that this will modify the Lithium abundances. 
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INTRODUCTION 

Many models of physics beyond the standard model 
have long-lived charged particles in their spectrum. For 
example, in the class of supergravity models where the 
gravitino is the lightest supersymmetric particle (LSP), 
the next lightest (NLSP) can be a long-lived charged slep- 
ton [l[ . Similar phenomenology arises in Universal Extra 
Dimension models [2]. Another example is the class of 
SUSY models with axino dark matter where again the 
NLSP could be long-lived [i, 0, I!- A third example 
is a supersymmetric model with almost degenerate LSP 
and NLSP such that the dominant decay channel for the 
NLSP is kinematically suppressed [1, 0| • 

Here we consider a new aspect of the early universe 
cosmology of these particles that affects the light element 
abundances. Other effects of these charged particles were 
OTeviously considered in [1, M, El, [l3, E, Q E H 

If the particle is negatively charged, it will eventually 
form a bound state with the positively charged nuclei. 
When the particle decays, the decay products will elec- 
tromagnetically interact with the nearby nucleus. If the 
nucleus is "'He, it may break up, creating ^He, T, D, p, 
n. These '^He, T and D nuclei will find other ^He nuclei 
and interact to form ^Li and ^Li . 

In addition to the production mechanism described 
above, charged particle decays result in a large non- 
thermal photon background that can destroy the created 
^Li and ^Li nuclei. When both the production and de- 
struction processes are taken into account, we find that 
the overall ^Li production could be in the same range 
as the Lithium isotopic abundance measured in few low 
metallicity popll stars [H, (H [H, [H, [H, [27 1 . 

Other effects can come from the elastic scattering of 
the charged particle produced in the decay from the 
"^He nucleus, and also from the change in the motion of 
the ^He in the bound state. We find that these effects 
produce insignificant changes to standard BBN abun- 
dances. 

In this work we only consider singly charged bound 
state of ''He nuclei [ii, Q [13, [IH . We ignore the pro- 
duction of neutral atoms of ''He bound to two X. We 
also do not consider the bound state of protons. This 



binding process starts below about 0.5 keV and could be 
important for lifetimes larger than about a year. Another 
interesting process we have not considered here is that of 
forming stable Be*-X bound states 28 1. A stable Be^-X 
bound state has important implications for the produc- 
tion of Carbon in the early universe. We hope to return 
to these issues in future work. 

While this work was being completed, two related arti- 
cles appeared on the arXiv. The three articles all consider 
different processes and hence complement each other. 
The first by Pospelov js^] pointed out that the nuclear 
processes with a photon in the final state will be mod- 
ified by the presence of the strong electric field of the 
X particle. Pospelov showed that the enhanced D on 
X-^He cross-section could produce orders of magnitude 
larger ®Li abundance. However, it is important to note 
that if the decay lifetime is about 10^ seconds or larger, 
then ^Li will also be destroyed due to the non-thermal 
photons from the decay. 

The second paper by Kohri and Takayama [s^l worked 
out the changes in the light element abundances due to 
the change in coulomb barrier in the presence of the , 
and due to the changes in the kinematics and energy 
balance for the reactions. They found that the Lithium 
abundances can be significantly modified from the stan- 
dard BBN predictions. 



RATES 

To calculate the magnitude of the effects we mentioned 
in the previous section, we need to calculate the following 
quantities. 

1. /{,, the fraction of Helium-4 nuclei bound to the 
charged particles. 

2. /p, the fraction of the bound ^He nuclei which are 
destroyed by photodisintegration. 

3. /a, the average fractions of '^He and T produced by 
the photodisintegration of a '^He nucleus. 

4. /i, the fraction of ■^Hc and T that formed in the 
photo-disintegration that find another *He nucleus 
and form ^Li . 
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FIG. 1: Solid curve is the fraction of helium nuclei that are 
bound to a long-lived massive charged particle. The number 
density of these charged particles is set by the requirement 
that the massive neutral particle it decays into be all the 
dark matter in the universe and the curve here is for a dark 
matter particle mass of 200 GeV. We have taken the decay 
lifetime to be much larger than the age of the universe at the 
temperatures shown above. 



In what follows we will concentrate on the 
^Li abundance. Some ^Li is also produced through 
'^He(^He, 7)^Be and ^He(T, 7)^Li. However, these reac- 
tions have an electromagnetic component and hence are 
slower than the ^Li production reactions. 

The total fraction of ®Li to Hydrogen (by number) is 
then hfphfLYp/A/{\ - Ypj^) where Yp is the total 
"'He mass fraction (bound or not) in the absence of the 
decay. The effects on the abundances of the other nuclei 
can also be similarly determined. 



CHARGED PARTICLE RECOMBINATION 

The bulk of ^He forms at T = r„ j - l^keV . If the 
decay happens significantly before this temperature j, 
then the effects we are interested in here are not impor- 
tant. For the rest of the discussion we will assume that 
the decay lifetime is larger than 1/ H{Taj). 

We also assume that the predominant component of X 
is negatively charged. Thus we are guaranteed that X 
will bind to the positively charged nuclei. 

The equation for the evolution of the fraction of 
charged X particles that are not bound (x/) is given by 



In the above equation, R is the recombination rate and 
(3 is the ionization rate. The recombination rate may be 
calculated in a manner analogous to the recombination of 
Hydrogen to yield R = 5{a/m)^/{BE/T) ln(BE/r). The 
ionization cross-section and hence /3 can then be found 
from the usual analysis of the equilibrium case. 



PHOTO-DISINTEGRATION OF HELIUM-4 

The photo-disintegration occurs through a one-photon 
exchange diagram, where the charged decay product 
emits a photon of invariant momentum = —Q^, which 
hits the Helium nucleus and breaks it up. 

Photo-disintegration can occur as long as the energy 
transferred to the hadronic system, is larger than the 
binding energy, which is here 28.8 MeV. Unfortunately, 
the cross-section for the process 7 + He^ any has not 
been measured for all the energies that we need. For large 
Q^, one can use the formulae for deep inelastic scattering, 
but these are only valid when is greater than Qq ^ 
IGeV^. 

To proceed, we will restrict our analysis to photon mo- 
menta satisfying > Qq- This will result in a conserva- 
tive estimate of the cross-section; the true cross-section 
will certainly receive additional contributions from pho- 
tons below this cutoff. 

For > Qq, the cross section is given by 
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Here ff{x) is the parton distribution function of flavor 
/ at longitudinal fraction x. Q/ is the corresponding 
charge. The final factor above is a kinematic constraint. 

Unlike the usual situation in deep inelastic scattering, 
here we should integrate over the allowed Q^. For a fixed 
X the range is Qq<Q^< xs. The allowed range of x is 
then between — and 1. 

s 

We then have to perform an integration over x. This 
is hampered by our lack of knowledge of ff{x). Quali- 
tatively, ff{x) increases at smaller x, so this integral is 
expected to be dominated by the smaller values of x. 

We will assume that e = <C 0.1, and restrict 
the range of the integral between x — e and x — 0.1. 
Over this ran ge, we can take approximately F2 {x) ~ 



^ff{x)Q)^im 



With these approximations, we find 
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Using IGeV-i - 
to about 27r x 10" 



' 0.2/m, the cross section evaluates 
4 X 0.04 X lQ-^^cm^ln{lQQl/s) ~ 
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0.3 X 10~^ln{lQQ'^/ s)'mb in almost exact agreement with 
a calculation from GEANT. 

To derive a probability from this quantity, we must 
multiply by the flux, and for this we must examine the 
geometry. The helium nucleus orbits around X at a 
distance given by d ~ mT^ jZa. = 3.6 fm where we 
have taken m to be the mass of the helium nucleus. 
The size of the helium nucleus may be taken to be 
r = 1.44^1/^ fm = 2.3 fm. The solid angle subtended 
by the nucleus at the position of X is 27r(l — cos (f)) where 
cos^ (j) = 1 — /(P and the probability of disassociation 



an 2(1 — cos cf) 
At:(P sin^ 6 cos^ ( 



0.012(T£, fm" 



(4) 



The probability of an interaction between the emitted 
particle and the nucleus is thus 0.012 x 3 x 10^^ = 3.6 x 
10~^ for a cross-section of 3 x lO^'^mb. 



NUCLEI FROM THE PHOTO-DISINTEGRATION 

The photo-disintegration of ^He described in the pre- 
vious section typically breaks up the nucleus. The nuclei 
will split into lighter elements "^He, T, D, p and n. The 
exclusive cross-sections for these processes are unknown. 
We will assume following [s^ that the helium nucleus 
is always destroyed, and produces '^He and T with the 
cross-sections 173 = ao = |crp. In other words, the pho- 
todisintegration produces a ^He or T nucleus in about 
half the interactions, and D in about half the interac- 
tions. We will assume that ^He and T are produced in 
equal amounts. 



PRODUCTION OF LITHIUM-6 

The ^He in the photo-dissociation of ''He are extremely 
relativistic. They can then scatter off background 
^He nuclei, and produce Lithium. The rate for this is 



given by r„ = 3.6 x lO-i^ypT^ ^y(a/40mb)(7?io/6)MeV 
where a is the appropriate nuclear cross-section. In com- 
parison, the expansion rate of the universe is H[T) 
2.5 X 10~'^'^T^^yMeV . Comparing the two rates, we see 
that the nuclear reaction rate will be larger than the ex- 
pansion rate of the universe until T ^ 0.3eV(40mb/ < 
av >). 

The nuclear reaction rate that we have quoted above 
assumes that the '^He are energetic enough to overcome 
the Coulomb barrier. This is certainly true for the nuclei 
at the moment of decay. However, these nuclei will lose 
energy due to a variety of processes [11], the most 
important of which are the energy losses due to Coulomb 
and Compton scattering. 

The time scale for these energy loss processes is much 
shorter than the nuclear reaction rate as well as the ex- 



pansion rate. As a conseguence one may write the prob- 
ability of ^Li to form as 



Eth 



dKa,+4K)v^nB{T) (5) 



where v is the relative speed, K is the kinetic energy 
of the nucleon, dK/dT is the total energy loss rate and 
(T3+4 is the cross-section for ^He on ^He reaction. Eth is 
the minimum energy required for the reaction to occur. 
The upper limit to the kinetic energy Kq is the energy 
with which the ^He nucleus is born. As it turns out, the 
integral is insensitive to the upper limit as long as it is 
larger than about 100 MeV. This is due to the fast energy 
degradation rate. 

This is not the end of the story. When the decay occurs 
the charged daughter particle carries off a large amount 
of energy. All of this gets converted to non-thermal pho- 
tons [36|. In the time it takes for these photons to ther- 
malize, they can find a ^Li (and other nuclei) and photo- 
dissociate it if the energy is sufficient to cross the thresh- 
old. 

We can calculate the non-thermal distribution of the 
photons following [1, [l^, [3^. The initial energy injection 
is rapidly brought to a quasi-static equilibrium state by 
pair production off of the background (thermal) photons 
and inverse Compton scattering. On a longer timescale 
(but still shorter than the age) processes like photon- 
photon scattering and Compton scattering redistribute 
the energy and lead to thermalization. We assume that 
the thermalization is mainly due to Compton scattering 
and calculate the differential non-thermal photon number 
density following [15 1. 

With the expression for the non-thermal photon distri- 
bution at hand, one may write the rate for ®Li production 
as 



dln(r) 



riBiT) 



tH{T) 

^6,7 



fL{T)hJph{T) 



H{T) 



(6) 



— where subscript "6" stands for ^Li and Fg.-y is the total 
photo-dissociation rate for 7-l-^Li anything. Y de- 
notes the number density relative to the number density 
of baryons. 

Similarly, the interaction of "^He and T with the back- 
ground protons and photons can degrade the nuclear en- 
ergy. We have checked that these are negligible effects 
for the temperatures of around 10 keV and lower that we 
are interested in. 

The non-thermal photons may also dissociate the 
■^He before they can collide with ''He nuclei. To estimate 
this, we note that the minimum threshold for '^He de- 
struction is about 5.5 MeV. In order to produce photons 
with this energy, one would have to be at a temperature 
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FIG. 2: Curve shows the ratio of number density of ®Li to 
that of H produced due to the process described in the text. 
The model assumed to make this plot is a gravitino Super- 
WIMP with mass of 100 GeV and a heavier charged slepton 
like stau. We start with an initial ®Li/H ratio of 1.4 x 10~^^. 



lower than about 2 keV. Thus, this could only be a prob- 
lem for long lifetimes. Also, the spectrum at these ener- 
gies is falling sharply 15|, |3g| . We estimate that this ef- 
fect cannot be significant unless decay lifetimes are larger 
than 10^ seconds and energy released is larger than TeV. 

The bottom-line from the above calculation is that it 
is possible to produce ^Li of order 10^^^ (relative to the 
number density of H) in large regions of decaying charged 
particle model parameter space. This is clear from the 
Figure [2] where we have assumed that a slepton is de- 
caying to a gravitino and lepton, with the gravitino hav- 
ing a mass of 100 GeV. As the mass of the gravitino is 
lowered, the amount of ^Li produced decreases. Note 
that the masses where a reasonable amount of ^Li is pro- 
duced corresponds to astrophysically interesting lifetimes 
of about a month [s^, . We note in closing that the 
dissociation by non-thermal photons will be insignificant 
for models where the mass of the decaying and daughter 
particles is fine-tuned to be small J5]. 



PRODUCTION OF LITHIUM-7 AND OTHER 
REACTIONS 

After the decay, all the "'He nuclei have the orbital ki- 
netic energy ~ 0.16 MeV. The kinetic energy is still not 
enough to initiate any endothermic nuclear reactions. We 
therefore consider the following processes '^He(a, 7)''Be 
and ^H(a,7)^Li. Proceeding as in the ^Li production 
case, we find that this process produces insignificant 



amounts of ^Li and ^Be . 

The ''He nuclei could also receive a large kick due to 
elastic scattering off of the decay produced charged parti- 
cle. Our calculations show that for small energy releases 
of order 10 GeV, it is possible to produce ''Li in quantities 
comparable to the standard BEN yield. However, in con- 
crete cases like the Super- WIMP scenario [l|, [l^ : the life- 
time for such small energy releases is large which means 
that the destruction rate due to non-thermal photons is 
large. The destruction effect dominates and significant 
amounts of ''Li are not produced. 

We also note that the break-up of ''He produces en- 
ergetic neutrons and protons that in turn can break up 
the background ''He nuclei. However, we have seen that 
only about 1 in a million bound '*He is broken up and 
hence the destruction due to these energetic protons and 
neutrons is negligible. 



DISCUSSION 

We have concentrated here on the collision of ^He with 
background ^He nuclei that would produce ^Li nuclei. 
Standard BBN does not produce much ^Li and thus the 
expectation is that this could be an important probe. 
The SBBN yield (without decays) is about 10~'^ com- 
pared to the number density of H. This should be com- 
pared to the measured lithium isotopic ratios in some low 
metallicity popll stars which (conservatively) range be- 
tween 0.01 and 0.1. A reasonable constraint on ^Li then 
is that its abundance should not be much larger than 
that of ^Li nuclei (about 10"^^° relative to the H number 
density). It is of course possible that the ^Li primordial 
abundance is larger, but then we are left with the ques- 
tion of why ^Li is depleted in much larger amounts than 
^Li by astrophysical mechanisms. 

We found that in some regions of parameter space it 
is possible to produce ^Li in quantities that match the 
observed "plateau" abundance. This does not imply that 
the decaying charged particle scenario can explain the 
^Li abundance. A full calculation including all the light 
elements would be necessary to ascertain that. We also 
note that there are astrophysical scenarios wherein the 
observed abundance of ^Li finds an explanation (e.g., (sol . 

3, EH, El H El El). 

The processes we have described here are clearly rele- 
vant for determining the light element abundances. Much 
more work needs to be done before these processes can 
be combined with previous analyses to put bounds on the 
model parameter space. Throughout this work we have 
neglected the bound states of other light elements. The 
fraction of bound states of D, T, '^He , ^Li and ^Li are 
small owing to the much smaller abundance of the respec- 
tive light elements. However, they could still have signif- 
icant effects because of enhanced nuclear cross-sections. 
The correct treatment of this problem including all the 
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cross-sections and bound states is beyond the scope of 
present work. 

In summary, we have studied an aspect of the decay of 
charged particles in the early universe. We looked at the 
electromagnetic bound states of the charged particle with 
^He nuclei and considered the effects of the decay on the 
light element abundances. We found that cosmologically 
relevant quantities of ^Li could be produced as a result 
of these decays. 
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